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In this paper we study the resolution ol images illuminated by sources composed of TV + 1 photons 
in which one non-degenerate photon is entangled with A'^ degenerate photons. The A'' degenerate 
photons illuminate an object and are collected by an A'^ photon detector. The signal from the 
photon detector is measured in coincidence with the non-degenerate photon giving rise to a ghost 
image. We discuss the case of three photons in various configurations and generalize to A'^-l- 1. Using 
the Rayleigh criterion, we find that the system may give an improvement in resolution by a factor 
of A" compared to using a classical source. For the case that the A'^-photon number detector is a 
point detector, a coherent image is obtained. If the A'-photon detector is a bucket detector, the 
image is incoherent. The visibility of the image in both cases is 1. In the opposite case in which 
the non-degenerate photon is scattered by the object, then, using an A'-photon point detector may 
reduce the Airy disk by a factor of TV. 

PACS numbers: 42.50.Dv, 42.30.Kq, 42.50.St, 07.07.Df 

I. INTRODUCTION 

Diffraction puts a limit on the the resolution of optical devices. According to the Rayleigh criterion [TJ [5], the 
ability to resolve two point sources is limited by the wavelength of the light. The Rayleigh or diffraction limit is not 
an absolute limit and proposals to exceed it have been known for a long time [2] . Recently, new proposals to improve 
resolution beyond the Rayleigh limit have been made based on the use of entangled sources and new measurement 
techniques. Improving the resolving power of optical systems beyond the diffraction limit not only is of interest to 
the fundamental research, but also holds promise applications in remote sensing and quantum sensors. 

Classical imaging can be thought of as a single photon process in the sense that the light detected is composed of 
photons each of which illuminates the object, consequently, the image can be constructed one photon at time. What 
we mean by referring to this as classical is that the source of the light may be described by a density matrix with a 
positive P-function jSJU]. In this sense the Rayleigh limit may be thought of as a single photon limit. Recall that ideal 
imaging is a process in which there is a point-to-point mapping of the object to a unique image plane. Diffraction 
causes each point of the object to be mapped onto a disk, the Airy disk, in the image plane. 

One of the new approaches to improving resolution is based on using non-classical light sources. Quantum ghost 
imaging [5J El H E [3 HO] is a process that uses two-photon entanglement. The unique features of this process are 
that entanglement allows only one photon to illuminate the object while the second photon does not. All the photons 
that illuminates the object are detected in a single (bucket) detector that does not resolve the image. The point 
detectors that detect the second photon must lie in a specific plane. This plane is called the image plane although 
there is no image in that plane; the image is formed in the correlation measurement of entangled photons. The image 
is constructed one pair at a time. The resolution of this system has recently been discussed [iTJ [12]. Losses in this 
system affect the counting rate but not the quality of the image. 

A second approach using non-classical source is based on entangled photon- number states [TH], e.g., NOON state. 
When the number of entangled photons exceeds two there are many possible imaging schemes that can be envisioned 
and so the analysis of these cases is still being carried out. This interferometric approach achieves a sub- wavelength 
spatial resolution by a factor N and requires an TV-photon absorption process. Another quantum source used to study 
imaging is to generate squeezed states [T3]. The image can be reconstructed through the homodyne detection [T5] . 
However, both of these techniques are severely limited by the loss of photons. 

A second class of approaches to improving resolution uses classical light sources. One method uses classical light 
with measurements based on correlations similar to ghost imaging and the Hanbury-Brown and Twiss experiment 
[inillZ]- This method has the advantage of being more robust with respect to losses [THlIin]- Another approach is to 
build an interferometric lithography with use of classical coherent state [201 EI], which has similar setup to the case 
using entangled photon-number states. 
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In this paper we will consider improving spatial resolution beyond the Rayleigh diffraction limit using quantum 
imaging with an entangled photon- number state |1,-/V). In our imaging scheme by sending the N degenerate photons 
to the object while keeping the non-degenerate photon and imaging lens in the laboratory, a factor of A'^ improvement 
can be achieved in spatial resolution enhancement compared to classical optics. The assumptions required for the 
enhancement by a factor of N are that the N photons sent to the object scatter off the same point and are detected by 
either an A^- photon number detector or a bucket detector. This sub-Rayleigh imaging resolution may have important 
applications in such as improving sensitivities of classical sensors and remote sensing. We emphasize that it is the 
quantum nature of the state that offers such sub-wavelength resolving power with high visibility. However, the system 
is very sensitive to loss. While we give general results, our main concern will be with the case in which the object is 
far from the source and the detectors and optics are close to the source. A different but related approach to the one 
discussed here is given in [^ . 

We organize the paper as follows. We will discuss our imaging scheme with entangled photon- number state |1,2) 
in some detail in Sec. II. In previous work [2^ I24j we have shown that imaging occurs in correlation measurement, as 
in the ghost imaging case. Here we will show that under certain stringent conditions, the resolution can be improved 
by a factor of 2 compared to classical optics. In Sec. HI we generalize the scheme to the |1, A^) case and show that 
resolution improvement by a factor of A'^ can be obtained. In Sec. IV some discussions will be addressed on other 
experimental configurations. Finally we will draw our conclusions in Sec. V. In an appendix we discuss the meaning 
of the approximation that the A^ photons illuminate the same point on the object. 



II. THREE-PHOTON OPTICS 



We start with three photons because this is the easiest case to investigate the various configurations. Throughout 
the paper we shall assume that the source of the three photons is a pure state and that the three-photon counting 
rate for three point detectors is give by 

Rcc^^J^ dhj^ dhj^ dh\^{l,2,i)\\ (1) 

where the three-photon amplitude is determined by matrix element between the vacuum state and the three-photon 
state 

vI/(l,2,3) = (0|i?i+)£;^+)4+V), (2) 

and 

E^j^\p3^Z]ii])= j di^j J d^ajEjfj{ujj)e-"^^*^gj{aj,ujj;pj,Zj)a{dj,ujj), (3) 

where Ej = y'?ia;j/2eo, dj is the transverse wave vector, and a{dj,u)j) is a photon annihilation operator at the output 
surface of the source, 

[a{d,iLj),a^d/,uj/)] = S{d - d/)S{uj - uj/). (4) 

The function Jj{Lo) is a narrow bandwidth filter function which is assumed to be peaked at ilj. The function gj is 
the Green's function [2 , 7^ that describes the propagation of each mode from the output surface of the source to the 
jth detector at the transverse coordinate pj , at the distance from the output surface of the crystal to the plane of the 
detector, Zj. is referred to as the three-photon amplitude (or three-photon wavefunction) . 

We start with the case in which the source produces three-photon entangled states with a pair of degenerate photons, 
that is ■0 ^ "01,2 

IV'1,2) = j dujiduj2 j d^aid^a25{2uji+ UJ2- ^)5{2di + d2)a\d2,U2)[a\di,ijJi)]'^\Q), (5) 

where is a constant, a;i^2 and di^2 £^re the frequencies and transverse wave vectors of the degenerate and non- 
degenerate photons, respectively. The i5-functions indicate that the source is assumed to produce three-photon states 
with perfect phase matching. We assume the paraxial approximation holds and that the temporal and transverse 
behavior of the waves factor. The frequency correlation determines the three-photon temporal properties. The 
transverse momentum correlation determines the spatial properties of entangled photons. It is this wave-vector 
correlation that we are going to concentrate on. As discussed in [231, several imaging schemes can be implemented 
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FIG. 1: (color online) Schematic of quantum imaging with a three-photon entangled state |1,2). di is the distance from the 
output surface of the source to the object. Li is the distance from the object to a 2-photon detector, Di. d2 is the distance 
from the output surface of the source to the imaging lens with focal length / and L2 is the length from the imaging lens to 
a single-photon detector D2, which scans coming signal photons in its transverse plane. "C.C." represents the joint-detection 
measurement. 



with this three-photon source. To demonstrate spatial resolution enhancement beyond the Rayleigh diffraction limit, 
consider the experimental setup shown in Fig. [T] It will be shown that for this configuration the spatial resolving 
power is improved by a factor of 2, provided the degenerate photons illuminate the same point on the object and are 
detected by a two photon detector. 

As depicted in Fig. [l] two degenerate photons with wavelength Ai are sent to a two-photon detector (Di) after 
illuminating an object, and the non-degenerate photon with wavelength A2 propagates to a single-photon detector 
(D2) after an imaging lens with focal length /. The three-photon amplitude ^ for detectors Di and D2, located at 
and (22, P2), now is 

^ ^ ^1,2 = {0\Ei+\p2, Z2,t2) Zl, il)]'|Vl,2), (6) 

Following the treatments in [51 [71 [53], we evaluate the Green's functions gi{di,uj2', Pi, zi) and g2{a2,u!2', P2, Z2) for 
the experimental setup of Fig. [T] assuming that the narrow bandwidth filter allows us to make the assumption that 
ujj — flj + Vj where \vj \ <C and 2f2i + U,2 — ^■ 

In the paraxial approximation it is convenient to write 

then 

where we replace by fij in = = 27r/Aj, z\—d\ + Li, and Z2 = '^2 + L2, respectively. In Eqs. ([8I and 

y, on 

^-1,2 = e'(2^^i^i+"^^2)$i_2, (10) 



Xi(ai,r!i;pi,zi) = e^'^^ J d^poA{po)e^^e"'^T^ (8) 

X2(a2,"2;/32,22) = 6 IdpiC ^ f e 1-2'^"'^ (9) 



([9]), A{po) is the aperture function of the object, and po and pi are two-dimensional vectors defined, respective 
the object and the imaging lens planes. With use of Eqs. ^ and ([5|, the three-photon amplitude ([6| becomes 



where Tj = tj — Zj /c and 

$1,2 = J dvidv25{2v^ + I/2)e*^'"^"^+"^"^Vl(r!l + V^f f2{^2 + V2)B^,2- (11) 



where 

Bi.2 - So / d' poA{po)e'^^ e-'^^ I d' i^MoV^^ I d'pie^^'^'T'e 



'2 /^/-/XJ^I— I J- / ^2 « 2 (t- j) i-[-P2Pl 



d^aie-'^'^'^'^^ + '^'>e-'^'<^P'-P"-P^\ (12) 

where we collect all the slowly varying quantities into the constant Bq. To proceed the discussion, in the following we 
will consider two different detection schemes. One uses a point two-photon detector for two degenerate photons after 
the object and the other has a two-photon bucket detector. 
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A. Point Two-Photon Detector Scheme 



In this detection scheme, a point two-photon detector is necessary to retrieve the information of degenerate photons 
scattered off the same point in the object. We therefore make the key assumption that the detector Di is only sensitive 
to the signals from the same point in the object, i.e., 5{po — p^) [The validity of this assumption is addressed in the 
Appendix]. With this assumption, Eq. (12 1 becomes 



Bi. 



Pie 



Vl,2 f)p 



aie 



(13) 



Completing the integration on the transverse mode di in Eq. ( 13 1 gives 



Bi.2 = Bo I dVo^^(Po)e'^^'^^°'''^+ ''i+(^^^/^i)''2 'e-'^^^^ 



„ ,- ^2 I Pj I r 1 1 ii .[Pi I Po 1 

d Pie ^ ll-a ^d2 + (^l/2A2)<il f'g '^^f LL2^d2 + (Ai/2A2)dil 



By imposing the Gaussian thin-lens imaging condition in Eq. (14) 

1 

d2 + {Xi/2X2)di' 
the transverse part of the three-photon amplitude reduces to 



1 

7 



1 

L2 



Bi, 



Bo 



dl + (2A2/>.i)d2 ' 



somb 



2ttR 



A2M2 + (Ai/2A2)di 



Po 



P2 



(14) 



(15) 



(16) 



where R is the radius of the imaging lens, R/[d2 + (Ai/2A2)(ii] may be thought of as the numerical aperture of the 
imaging system, and m = i2/['^2 + (Ai/2A2)di] is the magnification factor. In Eq. (16) the Airy disk is determined, 
as usual, by somb(a;) — 2Ji{x)/x, where Ji{x) is the first-order Bessel function. 

Before proceeding with the discussion of resolution, let us look at the physics behind Eqs. (15) and (16). Equation 



( 15 1 defines the image plane where the ideal the point-to-point mapping of the object plane occurs. The unique point- 



to-point correlation between the object and the imaging planes is the result of the transverse wavenumber correlation 
and the fact that we have assumed that the degenerate photons illuminate the same object point. Let us make a 
comparison with the two-photon and three-photon geometrical optics [71 111! 123] . In the Gauss thin lens equation 
the distance between the imaging lens and the object planes, d2 + {\i/2X2)di is similar to the form that appears in 
the non-degenerate two- photon case except for the factor of 2. This factor 2 comes from the degeneracy of the pair 
of photons that illuminate the object. As we will show below, this factor of 2 is the source of the improved spatial 
resolution. Equation (16) implies that a coherent and inverted image magnified by a factor of m is produced in the 
plane of I?2- Of course, there really is no such image and the true image is nonlocal. The point-spread function in 



Eq. ( 16 ) is generally determined by both wavelengths of the degenerate and non-degenerate photons. 



To examine the resolution using the Rayleigh criterion, we consider an object consisting of two point scatters, one 
located at the origin and the other at the point a in the object plane. 



MPof = Al6{po) + Al6{po - d). 



By substituting Eq. (17) into (16) we obtain 



Bi.2 = Bo A^somb 



2t:R 



P2^ 
L2 



s^'^M^somb 



2nR 



P2_ a 

L2 d2 + (Ai/2A2)di 



where the phase 



'P2 = Ki 



di + d2{2\2/Xi) 



a - {Pi + A) 



(17) 



(18) 



(19) 



indicates that the image is coherent. For a point 2-photon detector, we require pi = ff^ in Eq. (19 1. As is well-known 
[5] for coherent imaging the Rayleigh criterion is not the best choice for characterizing the resolution, however, it is 
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indicative of the resolution that can be attained and it is convenient. For a circular aperture, the radius of the Airy 
disk, ^, is determined by the point-spread function, which is 



C = 0.61- 



R 



(20) 



Note that the radius of the Airy disk is proportional to the wavelength of the non-degenerate photon. This is the 
standard result as obtained in classical optics. Using the Rayleigh criterion, the image of the second term in Eq. (18 1 
is taken to lie on the edge of the Airy disk of the first term, therefore. 



0.61- 



R 



2A2 



di 



(21) 



We see from Eq. (21) that the resolution depends on the wavelengths of the degenerate and the non-degenerate 
photons. In the case that di » ^2, so that ^2 + (Ai/2A2)di, is approximately (Ai/2A2)(ii. In this case Eq. (15) implies 
that L2 ^ f and the radius of the Airy disk approaches to 1.22X2/ / R, and 



= 0.61 



Aidi/2 
R ■ 



(22) 



Equation (22) shows a gain in spatial resolution of a factor of 2 compared to classical optics. Furthermore, there is 
no background term which is characteristic of the quantum case. 



B. Bucket Detector Scheme 



If the two-photon detector is replaced by a bucket detector and the two degenerate photons are collected by two 
single-photon detection events, located at {Li,pi) and (Li,//i), in the bucket, Eq. (12) becomes 



Bi.: 



(23) 



Under the assumption that the two degenerate photons are scattered ofi^ the same point in the object, Eq. (23) 



exp 



takes the similar form as Eq. (13), except that the second phase term in the first integrand of (13) is replaced by 



] . It is easy to show that the Gaussian thin- lens equation takes the same form as Eq. (15). By 



performing the same analysis as done in Sec. IIA on the resolving two spatially close point scatters, the three-photon 
amplitude ( 18 1 now is 



Bi,2 = Bo( Agsomb 



27ri? 



A, 



P2_ 

L2 



e^'^Mlsomb 



2ttR 



A, 



P2_ a 

L2 d2 + (Ai/2A2)di 



(24) 



Since the bucket detector gives no position information, we must square the amplitude and integrating over the bucket 
detector. 



/ - / d^p, / d'p[\B,,2\^ - sl\Bo\^ lAolSomb^ 



.f2TTR 



V A2 



P2_ 
L2 



27rR 



A, 



P2 



L2 d2 + (Ai/2A2)di 



(25) 



where Sb is the area of the bucket detector. It is easy to see that the spatial resolution improvement is the same as in 
Sec. IIA, the difference is that now we get an incoherent image. The advantage is that a two photon bucket detector 
should be easier to construct than a point two photon detector. 



III. TV + 1 PHOTON OPTICS 



In Sec. II, we have shown that with the entangled photon-number state |1, 2), the ability to resolve two point sources 
in the object can be improved by a factor of 2 by sending two degenerate photons to the object while keeping the 
non-degenerate photon and imaging lens in the laboratory. In this section, we are going to generalize the experimental 
configuration (Fig.fTl) with use of the entangled state of |1, iV), as described in Fig. [2] For simplicity, we first address the 
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FIG. 2: (color online) Generalization of quantum imaging with A*" + 1 entangled photons in state \1,N). For notations please 
refer to Fig. [T] except that here Di is an A'^-photon detector. The image is formed in the coincidence measurement and is not 
localized at either detector. 



case shown in Fig. |2] where the N degenerate photons traverse to the iV-photon detector, Di, after the object and the 
non-degenerate photon propagates to the single-photon detector, D2. The assumption required for the enhancement 
by a factor of N are that the N photons sent to the object scatter off the same point and are detected by the iV-photon 
detector, Di. 

The TV -|- 1 photons are assumed to be in a non-normahzed pure state 



|^/;ijv) = / duJiduj2 I (faid'^a26{NuJi+uj2-^)S{Nai+d2)ai (at )^|0). 



(26) 



Again the (^-functions in Eq. (26) indicate perfect phase matching. The N + 1-photon coincidence counting rate is 
defined as 



Rr.c — 



^ f dh f dt2--- f dtiV+l|«'l,7v(l,2, 
^ Jo Jo Jo 



(27) 



where ^'i^at is referred to as the N + l-photon amplitude. That is 

= {0\Ei+\p2,Z2,t2)[E[+\p^,z,,t,)]^\^/j^^N). (28) 

Following the procedure done for the |1,2) case, we calculate the transverse part of the N + 1-photon amplitude 
*i,Ar (p8| as 



N 



N 



B 



l.N 



- Bn 



d^PoA{po)e e" 



d'p'MV^^e-' — 



d Pie 



^2 S > Q 



N fold 
-iai-(Npi-po 



Po) 



(29) 



Here *&i,Ar(Ti, T2) describes the temporal behavior of entangled three photons. By applying the same argument that 
the iV- photon detector Di only receives the signals from the same spatial point in the object, Eq. (29) can be further 
simplified as 



N 



Bo I dVoA^(Po)e^^^e-""^ j d^pie^ 



K2\pir , 1 1 



/I 1 \ ■ -'^9 -' — 



(30) 



Performing the integration on the transverse mode a.\ in Eq. ( 30 ) gives 
Bx,N = Bo [ d^poA^{po)e 



^£il+(NA2/Ai)d2Jg -Li 



, „ j '<2\Pl\ r 1 I 1 11 _ 

X / d Pie ^ l-L2^<i2 + (Al/iVA2)di /Ig 



(31) 
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where, again, we have assumed multimode generation in the process. Applying the Gaussian thin-lens imaging 
condition 



1 

7 



1 



1 



d2 + (Ai/7VA2)di 



the transverse part of the N + 1-photon amplitude (|31| between detectors Di and D2 now becomes 

"27ri? 



1 1 - NKipi-po 

ll^ ' <ii+(N \2/ ^1)^2^ e ^1 somb 



A2 



P2_ 
L2 



Po 



d2 + {\ilN\2)di 



(32) 



(33) 



As expected, Eqs. (32l and (33l have the similar forms as Eqs. (151 and (I6I for the |1,2) case. The unique point-to- 
point relationship between the object and the imaging planes is enforced by the Gaussian thin-lens equation (32 1. The 
coherent and inverted image is demagnified by a factor of L2/[d2 + (ii(Ai/A^A2)]. The spatial resolution is determined 
by the width of the point-spread function in Eq. ( 33 1 . Note that a factor of N appears in the distance between the 
imaging lens and the object planes, d2 + di(Ai/A^A2). We emphasize again that the image is nonlocal and exists in 
the coincidence events. 

To study the spatial resolution, we again consider the object represented by Eq. (17 1. Plugging Eq. (17 1 into (33 1 
yields 



N 



Bo { A^somb 



/2nR 



For N single photon detectors located at p| 



P2 

L2 



e''^"Al^somb 



f^i^' the phase is given by 



2t:R 



P2_ 

L2 



{Xi/N\2)di 



(34) 



<Pn = Ki 



di+d2{N\2/Xl) 



«-(pf^+pr + ---) 



Li 



(35) 



For a point A^-photon number detector, we require /5| 



-(2) 
Pi 



and a coherent imaging is achievable in this case. 



The first term on the right-hand side in Eq. ( 18 1 gives the radius of the Airy disk, which is the same as the |1, 2) case, 
see Eq. (20 1. Applying the Rayleigh criterion, the minimum resolvable distance between two points in the transverse 



plane now is 



0.61- 



A2 
R 



NX. 



-di 



For the case of = 2, Eq. (36 1 reduces to Eq. (21 1. In the case that di ^ ^2, this becomes 

Xidi 



flni = 0.61 - 



NR 



(36) 



(37) 



As expected, Eq. (37) shows a gain in sub-Rayleigh resolution by a factor of N with respect to what one would obtain 
in classical optics. We therefore conclude that in the proposed imaging protocol, the spatial resolving power can be 
improved by a factor of A^ with use of the entangled photon- number state |1, TV). Furthermore, because we are using 
an entangled state with a specific type of detector, the image has high contrast because of the lack of background 
noise. 

By following the analysis in Sec. IIB, we can show that by replacing the A^-photon detector with an A^-photon 
bucket detector, we get an incoherent image but the sub-Rayleigh imaging process is not changed. 



IV. DISCUSSIONS AND OTHER CONFIGURATIONS 



In the previous two sections, we have analyzed a novel ghost imaging by sending A^ degenerate photons to the 
object while keeping the non-degenerate photon and imaging lens in the lab. We find that if the distance between 
the object plane and the output surface of the source is much greater than the distance between the imaging lens 
and the single-photon detector planes, we can gain spatial resolution improvement in the object by a factor of N 
compared to classical optics. In the cases that we have discussed in this paper, this enhancement beyond the Rayleigh 
criterion is due to the quantum nature of the entangled photon-number state. The assumptions required for such an 
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(a) 




(b) 



C.C. 




FIG. 3: (color online) Other schematics of quantum ghost imaging with three entangled photons in state |1,2). (a) Both the 
imaging lens and the object are inserted in the non-degenerate photon channel, (b) The imaging lens is placed in the degenerate 
photon pathway while the object is in the non-degenerate optical pathway. 



enhancement are that the N degenerate photons sent to the object scatter off the same point and are detected by 
either an TV-photon number detector or a bucket detector. An A^-photon bucket detector is much easier to realize 
than an iV-photon point detector. Such a bucket detector could be an array of single photon point detectors which 
only sent a signal to the coincidence circuit if exactly N of them fired. 

Besides the favorable configuration discussed above, one may wonder what happens if we switch the N degenerate 
photons to detector Di and the non-degenerate photon to D2 after an imaging lens and an object? Do we gain 
any spatial resolution improvement? To answer the questions, let us look at the |1,2) case as illustrated in Fig. |3ja). 
Following the treatments in Sec. IIA, after some algebra we find that the transverse part of the three-photon amplitude 
^ is 



Si,2 = Ba I d'poA{po)e 

,2 ' — 2 — [3- + 



2 ^ d'> 



2 e 



X 



d Pie 



d2+(Ai/2A2)ti /J 



-iK2pr 



d2+(Ai/2A2)ti 



(38) 



In the derivation of Eq. (38 1, the Green's functions associated with each beam give 

Xi(ai,f^i;pi,ii) = e~ 
X2{a2,^2;P2,Z2) = eT 



'1 e 



ipi-ai 



2 if 2 



K2IPLI1/ 1 1 



e 2 . 



d Pie 



Applying the Gaussian thin-lens imaging condition 



1 



1 _ 1 

d2 + (Ai/2A2)Li " 7' 



the transverse spatial part of the three-photon amplitude (p58| reduces to 
-Si, 2 = Bq d^poA{po)e 



^2 d'^i ^ I j,^ somb ( 



- ( 7^ + TTT ) ,■ -^'2P2-P. 

2 e 2 



V A2 



Po_ 

d' 



Pi 



d2 + (Ai/2A2)Li 



(39) 



(40) 



From this we see that the magnification is m = [^2 + (Ai/2A2)ii]/d2- Comparing Eqs. (39l and (40l with Eqs. ( [T5| 
and (16 1, we see that the distances between the object and the thin lens and between the thin lens and the imaging 
plane are interchanged. Since the degenerate photons are measured at the imaging plane in the setup of Fig. [3ja), 
the requirement of a point A^-photon detector cannot be relaxed. 
Computing the spatial resolution as in Sec. II we have 



Bi,2 — Bo 



Agsomb 



/27r 



R 



Pi 



d2 + (Ai/2A2)ii 



27ri? 



Pi 



d2 + (Ai/2A2)Li 



where ^' ^ K2[^-ii^ + jr) - ^] ■ The radius of the Airy disk is 



(41) 



(42) 
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If Li ^ (i2, C ~* "^fl^^ (^)' ^° ^'^^^ ^'^^ width of the point-spread function shrinks to one half its value compared to 
the classical cases. Applying the Rayleigh criterion to see the minimum resolvable distance between two point sources 
in the object. From the second term of Eq. (41 ) the minimum distance turns out to be 



0.61- 



2^2 



R 



(43) 



which only is a function of the wavelength of the non-degenerate photon; therefore, no spatial resolution improvement 
can be achieved compared to classical optics. 

Finally, we consider the configuration shown in Fig. |3jb) which was analyzed in [23 where it was shown that no 
well-defined images could be obtained. 

It is straightforward to generalize the above two configurations with use of the |1, N) state. By replacing the source 
state by the state |1, N) in Fig. [3]Ja), it can be shown that the radius of the Airy disk becomes 



0.61- 



A2 / Ai 



R ViVA. 



(44) 



"•^p^^ ij^), so the Airy disk shrinks to one iVth of its radius compared to classical optics. However, 



If il > d2, e 

if Li ^ d2, Eq. (44) gives the same result as in classical optics. Replacing the source with photon state |1,A^) in 
Fig. ^h), the above conclusion is still valid. The analysis has been presented in [2? and we will not repeat here. 



V. CONCLUSIONS 



In summary, we have proposed a quantum-imaging scheme to improve the spatial resolution in the object beyond 
the Rayleigh diffraction limit by using an entangled photon-number state |1,-/V). We have shown that by sending 
the N degenerate photons to the object, keeping the non-degenerate photon and imaging lens in the lab, and using 
a resolving A^-photon detector or a bucket detector, a factor of N can be achieved in spatial resolution enhancement 
using the Rayleigh criterion. The image is nonlocal and the quantum nature of the state leads to the sub-Rayleigh 
imaging resolution with high contrast. We have also shown that by sending the A'' degenerate photons freely to a point 
A-photon detector while propagating the non-degenerate photon through the imaging lens and the object, the Airy 
disk in the imaging can be shrunk by a factor of A^ under certain conditions. However, it may be possible to show that 
a similar effect can occur using non-entangled sources. In the language of quantum information, the non-degenerate 
photon may be thought of as an ancilla onto which the information about the object is transferred for measurement. 
Our imaging protocol may be of importance in many applications such as imaging, sensors, and telescopy. 
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APPENDIX A: VALIDITY OF THE ASSUMPTION MADE IN EQ. (13| 



In going from Eq. (12 1 to Eq. (13 1, we have made an assumption that requires the detector Di is only sensitive to 



the scattered photons from the same spatial point in the object. This allowed us to collapse the A^ integrations over 
the object into a single integral. In this Appendix, we give an example of how this assumption may be satisfied for 
multi-photon scattering off the target. Our example assumes that each point of the object transmits or scatters the 
light with a random phase which satisfies Gaussian statistics. The result is that the visibility decreases. 

We start with the case of 2 -I- 1 photons. From Eq. (12) the integration over the transverse vector cJi, which gives 



/ ■ -f^l IpoI^ r 2 I 1 1 . KlPl i Po f ■ -f^l IPo I r 2 I 1 1 

_-^lPl,2-Po 1/ -I \ ■ ^1 Po'P'o f ■ ^2 IP; 1^ r 1 I 1 J^-i _ - - \ P2_\ Po+p'q ] 

Xe~^ Ll g^fPtPo) g*2[di+(2A2/Ai)d2] / d ^ L ^2 + ( A 1 / 2 A2 ) di /^g 2 Pi" I + 2d2 + (A 1 / A2 ) d 1 J ^ (^^]^) 
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where pi j- is a point at which a photon is detected on the bucket detector, each point of the amplitude has a random 
phase associated with its transmission amplitude and, as usual, all the slowly varying terms have been grouped into 
Bq . Using the the Gaussian thin-lens imaging condition ( 15 1 gives 



Bi,2 = Ba I d^PoA{po)e 

■ ^1 Pl,2 P 



di-|-(2A2/Ai)d2-lg 



.■-f^llPo 



il-l T di+(2A2/Ai)d2 ' 



xe 



Ti g«'?'(Po)g*2[di+(2A2/^i)d2i sombl 





P2 


A2 


L2 



Po + Po 



(A2) 



Generalizing to the case of + 1 , using the Gaussian thin-lcns equation d32^ 



Bi. 



N 



Br 



d^Po,iA{pos)e e" 



- -f^lPl,! -Po,l 



'<i\Po,n\ 



) • • • / d'^Po,NA{po.N)e =^ 



■ ^lPl,NPo,N 



1 e 



y^gi<p(Po,N ) 2[d-^+{\2/N\i)d2\ somb 



P2 



P+ 



L2 d2 + iXi/NX2)di 



(A3) 



where p+ = jfJ2f=iPo,j- 

To compute the counting rate we first calculate the magnitude square of the amplitude averaged over the random 
phases. Starting with the N = 2 case and assuming that the ensemble average, (•••), over those phases satisfies 
Gaussian statistics so that 



(e.[0(Po)+0(p-'J-0(p;')-0(pT)l) = S{p, - f/^)6{f/, - pT) + S{Po - pT)5(p'o - P"'), 



(A4) 



We have assumed that the correlation length of the random phase is sufficiently small so that the Gaussian distribution 
can be approximated by delta functions. We find 



= \Bo\' d'po d'p'^A{pM{f^,)\'^omh 



V A2 



Po + P'o 



L2 2d2 + {\l/\2)di 



(A5) 



When we integrate over the bucket detector, the first term will be a constant while the second term will give us a 
delta function in po times the area of the bucket detector, Sh- Equation (A5l reduces to 



rfV.i / d'pi,2{\Bi,2?)^C+\Bo\^sl{^^) I d^po\A{p,)\'sonih' 



2ttR 



V A2 



P2 



Po 



(Ai/2A2)di 



(A6) 



First note that for the second term is similar to Eq. (25), the difference being the term in parenthesis which is the 
ratio of effect of diffraction to the area of the bucket detector, it is essentially the inverse of the Fresnel number. 
Computing the constant, C, is generally difficult and depends in detail on the geometry of the object, we can obtain 
an upper bound on C quite easily. 



\C\<sl\Bo\' 



d'po\Aipo)\' 



(A7) 



consequently, the visibility will be much less than for the ideal case discussed above. From Eq. (A6I the second term 
is proportional to LiXi which implies that as this product increases the visibility increases, however, recall for the 
case of sensors Li ~ di, so as this term increases the minimum resolvable distance also increases. 

The generalization to the case of A^ + 1 photons is straightforward. The ensemble phase average now becomes 



exp 



N 



N 



N 



(A8) 



Pn r=l 



where the A^ degenerate transmitted or reflected photons acquire random phases (j){po,j) and Pn is the set of permu- 
tations of the numbers (1, • • • , A). In Eq. ( A8 1 there are N\ terms. We can show that 



Bo\^ / d^po,i--- / d'^po,N\A{po,i) ■ ■ ■ A{po,N)\^somh- 



( 



2ttR 



V A2 



P2_ P+ 

L2 N[d2 + iXi/NX2)di 



e^*Tr2^„iPi,--(p°.--Po 



(A9) 
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When we integrate over the bucket detector, we get a complicated result. Two terms are simple, the identity per- 
mutation gives a constant and the single cycle subgroup give an incoherent image with a resolution that depends on 
\\/N. These are the only terms for N — 2. The remaining terms will lead to terms which are essentially constant. 
For = 3 we get 



d^PM / dVi,2 / d'piM\Bi,3\')=C + 3\Bo\'sl{ ^ 



xsomb 
xsomb 



2t:R 



2/27ri? 



V A2 



P2_ 
L2 



P+ 



3[d2 + (Ai/3A2)di 



d2 + {\i/NX2)di 



\A{po,i)\' 



(AlO) 



From Eq. (AlO) the second term shows explicitly how the general terms will lead to a complicated average over the 
illuminated area of the object. This result shows that the image will have very poor visibility for large N, it is not 
certain whether there might be arrangement of detectors for the TV photons which will give better results. 
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